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Abstract 



We analyse vacuum gravitational "soliton" solutions with real poles in the cosmo- 
logical context. It is well known that these solutions contain singularities on certain 
null hypersurfaces. Using a Kasner seed solution, we demonstrate that these may 
contain thin sheets of null matter or may be simple coordinate singularities, and we 
describe a number of possible extensions through them. 



OO 
ON 

^ 1 Introduction 

tjT)! The inverse scattering (BZ) technique of Belinskii and Zakharov |I| is now well known. 

It is essentially a solution-generating procedure for producing exact vacuum solutions of 
Einstein's equations for space-times admitting two isometries. Starting from some initial 
"seed" solution, the technique is based on the construction of a "dressing" matrix which 
is a meromorphic function of a complex spectral parameter A. 

For the case in which a vacuum space-time admits two hypersurface-orthogonal space- 
like Killing vectors, the "gravitational solitons" corresponding to particular poles of the 
dressing matrix generally describe perturbations of the gravitational field which propagate 
like finite gravitational waves on some background. Here we consider gravitational soliton 
solutions in such space-times which correspond to real poles in the scattering matrix. 
As originally pointed out by Belinskii and Zakharov these solutions exist in regions 
which are bounded by null hypersurfaces on which singularities occur. It is the purpose 
of this paper to reconsider the character of these singularities and the possible extensions 
through them. 

Carr and Verdaguer |2| have considered soliton solutions in a Kasner background and 
have interpreted the solutions with real poles as inhomogeneous cosmologies with shock 
waves in which the solitons propagate away revealing the Kasner background. However, 
as shown by Gleiser [[J and Curir, Francaviglia and Verdaguer these solutions must 
contain thin sheets of null matter separating the various regions. Gleiser || has also 
described alternative matter-free extensions, while Curir, Francaviglia and Verdaguer || 
have considered a real pole of arbitrary degeneracy for the diagonal "soliton" solution of 
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Carmeli and Charach || which was shown to correspond to a real (degenerate) pole by 
Feinstein and Charach 0. 

Others 0-[lO have considered soliton solutions with real poles in which the seed 



metric is a nondiagonal vacuum Bianchi II space-time. In this case, it is similarly possible 
to remove the coordinate singularity on the null hypersurface, although there may again 
be an impulsive gravitational wave and a 5- function in the Ricci tensor on the shock front. 

In this paper, we reconsider the physical interpretation of some soliton solutions with 
real poles. As summarised above, these can represent gravitational shock waves in some 
cosmological background. However, the extension to the background is not unique even 
if the possibility of thin sheets of null matter is excluded. In particular, we consider the 
one-soliton solution with a vacuum Bianchi I seed. By extending this through the shock 
front in various ways, we construct a number of different global solutions. 



2 Real pole solitons with a Kasner seed 

According to the BZ technique [0, we consider a vacuum space-time with two hyper- 
surface-orthogonal Killing vectors. In the case in which the isometries are spacelike, the 
metric can be written in the form 

ds 2 = 2e' M du dv - gijdx^ (1) 

where the 2-metric g (= gij) is a function of the two null coordinates u and v, and 
has the determinant |g| = a 2 . Einstein's vacuum field equations require that a satisfies 
the 2-dimensional wave equation and so can be written in the form a = f(u) + g(v), 
where f(u) and g(v) are arbitrary functions. It is also convenient to introduce another 
function fl(u,v) which is harmonically conjugate to a and given by (3 = f(u) — giy). In 
this and following sections, we will assume that the two null coordinates u and v are 
future-pointing. 

Here, we start with an initial Kasner seed solution (denoted by a suffix zero) which 
can be written in the form (1) with 

\ _ Mq a u a v 



where p is an arbitrary parameter. This reduces to a form of the Minkowski metric when 
p = ±1. According to the BZ technique, we work with the matrix 

*° (M ' M)= i (a 2 + 2/3A + A 2 )(^)/ 2 J 

which satisfies the appropriate equation and the condition ^ (u,v,0) = go- 

We also restrict attention here to the case in which there is a single real pole given by 
A = /i, where 

/i = uj - (3 + \J(uj - (3) 2 - a 2 

and io is an arbitrary real constant. (For a single real pole, we note that the alternative 
expression fj, — lu — (3 — ^ {uj — (3) 2 — a 2 simply corresponds to a rotation of coordinates 
x — > y, y — > —x.) Clearly the single soliton solution is only admissible in regions of 
space-time for which (uj — (3) 2 > a 2 . (If there are more than one such regions, these will 
normally be disjoint.) These regions will be bounded by an initial cosmological curvature 
singularity which occurs when a = 0. We will therefore take a > 0. They will also be 



bounded by null hypersurfaces on which (u — (3) 2 = a 2 . According to common usage, we 
will refer to these as "shock fronts" . 

Since a 2 + 2f3X + A 2 = 2u\ we obtain that on the pole trajectory A = \i 

/(2 WA *)-a+p)/2 o \ 

*o(w)=( Q (2 ^ } _ (1 _ P)/2 ) • 

The procedure is now well known and, for the nondiagonal case in the region u — (3 > a, 
the new solution can be expressed as 

\a otfj, J 

where, after introducing a new function s(a, 0) such that e s = fi/a and the constant 
c = plog(2uj), the matrix P is given by 



2cosh(ps + c) V a~ p e ps+ \ 

It may be noted that a change in the values of the free soliton parameters corresponds to 
a rotation of coordinates. For the single soliton solution, this is not physically significant. 
However, when attached to other regions, changes in these parameters can represent 
soliton waves with different polarization. With the above expressions, the new solution is 
given by 



a ( Oi v cosh[(l + p)s + c] — sinhs 



coshQos + c) \ —sinhs a p cosh[(l — p)s — c] 



(3) 



3 Extending the solution to a Kasner background 

As noted, the above solution using (3) is only defined in the region of space-time in which 
u — (3 > a > 0. This region is bounded by the null hypersurface (3 + a = uo on which 
/j, = a. It can be seen that, on this boundary, s = and the 2-metric is the same as 
that of the seed. Thus, it would seem to be possible to join this metric continuously with 
the seed Kasner solution in the region a > uj — j3, thus forming a gravitational shock 
front as described by Carr and Verdaguer 0. It is also possible to include a one-soliton 
solution in the region in which f3 — uj > a > 0. The standard global interpretation of 
this solution is then as two disjointed one-soliton regions moving apart leaving an exact 
Kasner background. This forms a composite space-time as illustrated in figure 1 in which 
there is an initial curvature singularity when a = 0. 

However, although the 2-metric is continuous in this construction, it may have discon- 
tinuous first derivatives which could introduce an impulsive gravitational wave component 
on the shock front. 

More seriously, we note that, using (3), the new expression for M in the region uo — (3> 
a > is given by 

-m = C V® coshps + c) e _ Mo 
\Juj — a — (3yjui + a — (3 

where C is an arbitrary constant. This clearly introduces a singularity on the shock front 
(3 + a = uj. However this can be seen to be a coordinate singularity which can be removed 
by a particular choice of the functions f(u). For example, using (2), the singularity can 
be removed by the choice / = |(a; — u 2 ). 
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Figure 1: This represents a one-soliton solution, using a and (3 as coordinates, in which 
the solitons propagate away to reveal a Kasner background. The space-time is a combi- 
nation of two one-soliton regions (denoted by I a and K) each with a real pole uji = —z±, 
together with a Kasner region (denoted by K). 



It may be recalled that the BZ method is based on the assumption that a = f + g is 
the same for the soliton solution as for the seed metric that was used to generate it. Thus, 
with the choice / = \{uj — u 2 ), the seed that is used in this case is a particular form of 
the Kasner metric which contains a coordinate singularity on the null hypersurface u = 0. 
i.e. in order to generate a non-singular one-soliton solution directly, it is necessary to 
start with a seed solution in a form containing a coordinate singularity. 

In the above construction, the space-time has been extended through the null hyper- 
surface u = 0, on which \x = a, to a "background" region which is part of a Kasner 
space-time with the same parameter p as the seed solution. However, to avoid a coordi- 
nate singularity in this region, it is necessary to put a = |(it + v), which is equivalent to 
the choice / = |(w + uj) and g = \{y — u). In order for this to be continuous with a 
nonsingular soliton region, it is necessary for the soliton region to have / = ~(a; — u 2 ) 
and g = ~(v — w) (with the region I?, the same but with u and v interchanged). From 
this, it follows that a in the soliton regions (and the associated seed) is different to that 
in any extended "background" region, and that the background only has the same form 
as the seed solution after a coordinate transformation. 

In the soliton region l a we now have a — \{v — u 2 ). This has been made continuous 
with a Kasner background in which a = ^(u + v) across u — 0. However, it may be 
recalled that discontinuities in the derivatives of a across a hypersurface induce nonzero 
components in the Ricci tensor, and hence in the energy-momentum tensor. In this case, 
these are given by 

e M e M 
T uu = -- — [a u ]6(u) or T vv = -- — [a v ]6(v). 
owa ona 

It can thus be seen that the discontinuity in the derivative of a across u = in the 
above expressions gives rise to an impulsive component in the energy-momentum tensor 
corresponding to a thin sheet of null matter located on this hypersurface. This has been 
described explicitly elsewhere 0, ||]. Moreover, with this time orientation, the matter 
has negative energy density JED . 



We now have a one-soliton solution (3) with a single real pole lui = —Z\ in the region 
represented as l a in figure 1. This has been extended to a Kasner region (represented by 
K in figure 1) across the shock front u = —Zi, and another one-soliton region (represented 
by Ift) in which uj — (3 < —a has been added together with a second shock front v — Z\. 
The metric function a = f(u) + g(v) is determined by the following expressions in each 
region: 
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Figure 2: This represents the one-soliton solution in which the solitons propagate away 
to reveal a Kasner background in the u, v coordinates. An initial "cosmological" curvature 
singularity occurs when a = 0. 



K: u > —zi, v > zi 





la- U < -Zi, V > Zi 



f = - {u + zif} 9 = \v 



lb- U > —Zi, V < Z\ 




g = \[zi- {zi-vf} 



Clearly there are discontinuities in the derivatives of a across the junctions u = —z\ and 
v = Zi. Thus, as described above, this configuration must contain thin sheets of null 
matter on these hypersurfaces. 

It may also be noted that the curvature singularity at a = is located, in these 
coordinates, on v — z\ — (u + zi) 2 = for v > Z\ and on u + Z\ — {z\ — v) 2 = for v < z\. 
Thus, although this singularity is spacelike, it approaches a null limit at the junction 
point v = —u = Zi. This is illustrated in figure 2. 

It is appropriate to work with a null tetrad such that ii = e~ M / 2 u t i and rij = e~ M ^ 2 v 
Since M is continuous across the shock fronts in these coordinates, these null vectors are 
well behaved throughout the space-time. Using these, it can then be shown that (provided 
p 7^ 0, ±1) the Weyl tensor has non-zero components ^o, ^2 and \l/4 in both the soliton 
and Kasner regions. Moreover, these components are bounded near the shock fronts, 
although they are not necessarily continuous across the shock which may also have step 
changes in the polarization into the soliton regions (in the nondiagonal case). In addition, 
the shock fronts themselves may also contain impulsive components. 

4 Soliton solutions with distinct real poles 

In this section we re-consider the n-soliton solutions with distinct real poles which have 
been described briefly elsewhere [0 and interpreted as solitons moving apart leaving an 
exact Kasner background. To be specific, we concentrate on the case of a soliton solution 
with two distinct real poles uoi = —Z\ and U2 = —Z2, where Z2 > Z\. This can again be 
interpreted as a composite space-time having one- and two-soliton regions and a Kasner 
background as illustrated in figure 3 using a, [3 coordinates. In this case, the shock fronts 
can be taken to be u = —Zx, u = —Z2, v = z\ and v = Z2- 

As described in the previous section, singularities can occur in the metric coefficient 
e . For the two-soliton solution this is given by 




(a 2 - /ii)(a 2 - /ii/i2) 2 (a : 



2 -/ii)' 



g = 



Z\ z 2 

Figure 3: This represents soliton solutions with two real poles with a and (3 as coordi- 
nates. The solitons propagate away to reveal a Kasner background. The space-time is a 
combination of one- and two-soliton regions (denoted by I a , I;,, II a , II^,, II C ) with real poles 
uj\ = —z\ and 102 = —Z2, together with a Kasner region (denoted by K). 



An explicit transformation which removes each singularity individually has been given by 
Diaz and Gleiser [TTJ (for the general case of n distinct real poles). However, this can only 
be applied to one singularity at a time as the resulting metric becomes discontinuous on 
the next singular wavefront. Below, we introduce a different gauge in which the metric is 
continuous throughout the space-time. 

The two-soliton solution is appropriate in the regions II a , Ilf, and II C indicated in 
figure 2. It is then possible to extend the solution into regions described by the one- 
soliton solution (3) and then further to the Kasner solution (2) as illustrated. We can 
adopt two future-pointing null coordinates u and v that are defined globally. It is then 
possible to introduce a gauge such that f{u) and g(v), and hence a(u, v) and f3(u, v) are 
continuous everywhere. One such gauge is as follows: 



K: 



u > —Z\, 



V > z 2 , 



la- ~Zi >U> -Z 2 , V > Z 2 , 



U > —Zi, Z2 > V > Zl, 



Ha- -Z2 > U, 



V > z 2 , 



ll b : —z\ >u> -z 2 , z 2 >v>zi, 



II C : u > -z\. 



Zi > V, 



f = 










9 = 


¥ 








f = 








{u + zxf - k{u + z x f 


9 = 


¥ 








f = 










9 = 




- 




- v ) 2 + k(z 2 - v ) 3 ] 


f = 


§[- 


Z 2 




{u + z 2 )% 


9 = 


& 








f = 


u- 


Zl 




(u + z x f - k{u + z x f 


9 = 






>2 


- v) 2 + k(z 2 - v) 3 ] 


f = 










9 = 






>1 





where k = (z 2 — Zi)^ 1 — (z 2 — Zi)~ 2 . 

With these expressions, it can be seen that the metric is C° everywhere and that all 
coordinate singularities have been removed. However, it can also be seen that there are 
discontinuities in the derivatives of a on the null hypersurfaces u = —Z\, u = —Z2, v = Z\ 



Figure 4: This represents soliton solutions with two real poles in the u, v coordinates. 
An initial "cosmological" curvature singularity occurs when a — 0. 



and v = Z2 which therefore must contain thin sheets of null matter. The structure in u, v 
coordinates is illustrated in figure 4. 

It may be noted that the space-times in the soliton regions may be diagonal or non- 
diagonal. However, some soliton parameters must be continued across the shock fronts. 
To be explicit, the one-soliton regions must contain solitons whose parameters are con- 
tinued from the two-soliton regions II a or II C . The region lit,, however, must contain a 
continuation of the two solitons which extend back from regions l a and lb and therefore 
its parameters are predetermined from those of II a and II C . 

Clearly, this approach can be generalised to include an arbitrary number of distinct 
real poles. For solutions with more than one real pole, these solutions necessarily contain 
thin sheets of null matter with negative energy density on the boundaries between the 
various regions. 

5 Possible extensions without sheets of null matter 

We now consider the possible extensions of the one-soliton solution which do not involve 
thin sheets of null matter. In this case, no generality is lost in making a coordinate shift 
to put oj = so that the shock front of the soliton region is then given by u = 0. We also 
continue to use two future-pointing null coordinates u and v. 

We may initially consider whether or not it is possible to construct an exact one- 
soliton solution which has the same global structure as that indicated in Figure 1, but 
without the presence of thin sheets of null matter. For this, the regions I a and h would 
be essentially the same, but the extension would not be to a Kasner background. If, 
in the one-soliton region I a we put a = \{v 2 — u 2 ), it may be possible in the extended 
region to choose a gauge such that a = |(w 2 + v 2 ). However, it is shown in the appendix 
that such an extension is not possible. We therefore look for alternative extensions to 
the region I a . Since the shock front is null, any extension will be non-unique. In fact, a 
number of possibilities are readily available. 

It may first be noted that the one-soliton region is algebraically general. There will 
therefore be a gravitational wave component propagating towards and through the shock 
front. In the soliton region, there will also be a gravitational wave component propagating 
in the opposite direction parallel to the shock front. However, this component need not 
occur in the extended region and, in this case, the extension will be to a plane wave region. 
Other possibilities include that in which the extension is to another distinct soliton region, 
or to that containing an arbitrary gravitational wave component propagating parallel to 
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Figure 5: A one-soliton region can be extended to a plane wave region. Gravitational 
wave components originating in the singularity propagate into the plane wave region but, 
in that region, no wave components propagate along the lines u = const. A non-scalar 
curvature singularity occurs on v — 0. 



the shock front (although only the linear case will be constructed below). 
5.1 A plane wave extension 

The simplest extension across the shock front u = is to a plane wave region in which 
/ = and g and M are functions of v only. The general structure is illustrated in figure 5. 

For the moment, let us continue to adopt a gauge such that g(v) = \v. In the extension 
to a plane wave region, we may now set /(it) =0 and a = —[3 = \v and the metric is 
given by 



e~ M = cM p2 - 1)/2 , g = - 








2 V (§)-*>, 



It is then appropriate to make the coordinate transformation 

v = v 2 /^ 

together with a rescaling to remove an unwanted constant so that the line element becomes 

ds 2 = 2du dv - £ 2 /(P 2+1 ) (>/( P 2 +l) dx 2 + ~-2 P /( P * + l) dy 2^ 

We can then make the further coordinate transformation 

l+p 

X = v 1+ p' 2 x 
-p 



Y = v^+p' 2 y 

, . s l+2p-p 2 , , , l-2p-p 2 



to put the line element in the form 



ds 2 = 2dv dr - v-\X 2 - Y 2 )dv 2 - dX 2 - dY 2 



which is the familiar form of a plane wave with amplitude profile h(v) = jj^jrp v 2 ■ Thus, 
except for the cases in which p — 0, ±1 in which the extension is to a flat region, the plane 
gravitational wave amplitude in the extended region is clearly unbounded when v = (or 
v — 0). This null hypersurface may then reasonably be considered to form a boundary of 
the space-time extended from the one-soliton solution as described in figure 5. 

It may be observed that this situation is qualitatively identical to the time reverse 
of a colliding plane wave space-time in which a future curvature singularity is formed 
following the interaction of initially plane waves. As in that case, the non-scalar curvature 
singularity that occurs in the plane wave region on v — can be interpreted as a "fold 
singularity" as described by Matzner and Tipler []T2| (see also 



s 



5.2 A soliton extension 



An alternative extension can be achieved by matching possibly different one-soliton solu- 
tions on either side of the shock front. On either side we can adopt the same gauge with 
/ = —\u 2 and g = ^v 2 , so that 

a = \( v 2 -v 2 ) and = -\(u 2 + v 2 ). (4) 

This clearly satisfies the required inequality for the one-soliton solution. Taking u to be 
time-oriented, we now have two regions. An initial region with u < and a second region 
with u > 0. 

It is also appropriate here to introduce alternative coordinates t and p where 

t= 72^ u + v ) and P = 72( V ~ U )' 

so that a = tp and p = t 2 . 

The curvature singularity at a = now occurs both when t = and when p = so 
that the space-time is defined only in regions for which v > \u\. However, this includes 
two regions in which the conditions for a one-soliton solution are satisfied. These occur 
on either side of the shock front as illustrated in figure 6. Of course, there is no reason for 
solitons in these two regions to have the same parameters, and it is possible to construct 
a compound space-time composed of two different one-soliton solutions joined across the 
shock front. For example, we can choose the soliton in one region to be diagonal, and that 
in the other non-diagonal. In addition, since / and g are taken to have the same form 
in both regions, there is no discontinuity in the derivatives of a, and so there will be no 
sheets of null matter across the shock front. This matter-free extension of the one-soliton 
solution was initially pointed out by Gleiser || at least for the case when p = —1. 

As a particular example, let us consider the case of the one-soliton solutions that are 
generated from the plane symmetric (type D) Kasner solution for which p = 0. The 
non-diagonal case is given by 

/ cosh s — sinh s \ 
g = tp\ where e = t/p. (5) 

\ — sinh s cosh s J 

This solution can be adopted in the region t > p > 0. Choosing C = —uj, the remaining 
part of the metric is given by 

2e- M dudv = dt 2 -dp 2 . (6) 

Let us now attach this to a background region u < and v > (i.e. p > t > 0) which 
is composed of the diagonal one-soliton solution with ap = Kasner seed, given by 



t/p 



This can be expressed in the form 

ds 2 = dt 2 _ df) 2 _ ^2 _ p 2 d0 2 (?) 

which after the coordinate transformation 

T = t cosh z, X = p cos 0, Y = p sin 0, Z = t sinh z 



Figure 6: This space-time is composed of two (possibly different) one-soliton solutions 
in the two regions denoted by I a , K. It is defined only for v 2 — u 2 > 0. Generally, there 
are curvature singularities on the hypersurfaces t — and p = 0. 



is clearly that part T > \Z\ of the Minkowski space-time 

ds 2 = dT 2 - dX 2 - dY 2 - dZ 2 

in cartesian coordinates. In this case, it can be seen that the junction between the two 
regions, which is the null hypersurface u = (or t — p), is an expanding sphere given by 
T 2 — X 2 — Y 2 — Z 2 = 0. The solution therefore describes a gravitational wave with an 
exact spherical wavefront propagating into a Minkowski background. 

It may also be observed that the metric (5-6) may be written in the form 

ds 2 = dt 2 - dp 2 - \t 2 (dz - d<pf - \p 2 (dz + d(j)) 2 

which is clearly simply a rotation of the flat metric (7). Thus, the space-time inside 
the spherical wavefront is just another version of Minkowski space. The single soliton 
solution obtained from a plane symmetric Kasner seed is a flat space-time. The only 
possible nonzero components of the Weyl tensor arise from the discontinuities in the 
derivatives of the metric across the wavefront. The solution described above therefore 
represents an impulsive spherical wave propagating in a Minkowski background. This is 
in fact a special case of the impulsive spherical gravitational wave that was constructed 
by Penrose |T3j using a "cut and paste" method. 

It may also be observed in this case that t = simply corresponds to a coordinate 
singularity in the Minkowski background. It is thus possible to add the Minkowski region 
in which < T < \Z\, and then to add the time-reverse of the solution for T < and 
t < 0. The global solution then describes a contracting impulsive gravitational wave in 
a Minkowski background which collapses to a point. After the collapse, the wave then 
re-expands as an exact spherical impulsive wave. The only singularity occurs at the event 
at which the spherical wave has zero radius. 



It may also be pointed out that Gleiser, Garate and Nicassio \T5\ have obtained a 
similar singularity-free solution in a one-soliton solution generated from a Bianchi VIo 
seed. In that work they describe the soliton perturbation as "erasing" the "cosmological" 
singularity that occurs in the seed. 

For the case of a more general Kasner seed, solutions can be constructed as outlined 
above. For p ^ 0, they will have curvature singularities both when t = and when p = 0, 
so that the complete space-time is the region t > 0, p > 0. However, these cases do not 
strictly have an axis at p = and the space is not asymptotically flat. 

5.3 A non-soliton extension 

In the plane wave extension of §5.1, the gravitational waves propagating in the direction 
v = const, simply continue into the plane wave region in which no waves propagate 



in 



in the opposite direction. However, it is natural to consider the possible existence of 
gravitational waves propagating in the extended region parallel to the shock front. Such 
components occur in the soliton extension of §5.2 in which both regions are algebraically 
general, although only specific wave profiles are permitted. We now consider an alternative 
extension which includes arbitrary gravitational wave components in the extended region 
- at least for the diagonal case in which such waves have constant aligned polarization. 
Such an extension can be constructed as follows. 

We start by choosing coordinates in the one-soliton region with / = — \u 2 and 
g = | v 2 so that a and j3 are given by (4). We then adopt the same expression for a(u, v) 
in the extended region so that there are no Ricci tensor components across the shock 
front. The metric in the extended region can then be taken in the form 

ds 2 = 2e- M dudv - a{e v dx 2 + e~ v dy 2 ), 

where 

V(u, v ) = p log a + V(u, v) 

and V is an arbitrary function satisfying V"(0, v) = 0, so that the metric remains contin- 
uous across the shock front. 

An expression for V satisfying the above property is given by the Rosen pulse solution 
| TBI which, in this context, takes the form 



V 



F(a) da f° F(a) da 



o yja - fy/a + g J-u 2 y/a + u 2 y/a + 



where F(a) is an arbitrary function. In the present situation we require F(Q) = and the 
continuity of F at a = must be determined very carefully. However, since a complete 
solution cannot be determined for the solution in this form, its application here is very 
limited. 

On the other hand, an explicit representation for V for which a complete solution can 



be obtained, and which satisfies the required properties above, has been given in ||17||-||19||. 
These solutions are expressed as a sum over explicit components, each of which have the 
self-similar form 

V h (f,9) = (f + 9) h Hk('££ 



J+g, 

where k is an arbitrary real parameter. Putting ( = —/3/a, the functions iifc(C) satisfy 
the linear equation 

(C 2 - l)H'l + (1 - 2k)CH' k + k 2 H k = 0, 
together with the initial condition Hk(l) = 0, so that they satisfy the recursion relations 

H k (C) = £ if fc _i(0 d£ or H' k (C) = H k ^(C). 

These solutions can be expressed in terms of standard hypergeometric functions in the 
form 



where, for integer k, c k = (— l) /c 2 A T(|)/r(fc + |). The papers however, are con- 

cerned with gravitational waves with distinct wavefronts propagating into certain simple 
backgrounds. For such situations, it is only necessary to consider cases in which k > \. 
However, in the present situation, it is necessary to choose the minimum value of k in 
order to remove the singularity that would otherwise arise in the metric coefficient e _M , 
and it can be seen that this must involve the component k = 0. 



I I 



In determining an explicit solution for the extension, we may use the fact that 



H (C) = cosh" 1 C and #-i(C) 



It is convenient initially to consider the solution V = aoHo((), where ao is a constant, so 
that 

V = ploga + a Ho(() 

= p\og(f + g) + a H (^). 

The remaining vacuum field equations as given in [|17]] can then be integrated, yielding 

(£2 _ l)a*/2 

= C\f' 9 '\(^ + V^) 2ao \f + g)^ 2 -^ 

where C is an arbitrary constant. With the above expressions for / and g, it can imme- 
diately be seen that a coordinate singularity is avoided only if a — ±1. Taking a = 1 
and a particular value for C, we obtain that 



v — uj 

indicating that the metric is now continuous as required across the shock front. However, 
it may be noticed that this case in which V = -ffo(C) i s J us t the diagonal case of the one- 
soliton solution, and therefore belongs to the class of extensions discussed in the previous 
subsection. 

A more general extension in which the metric is diagonal in the extended region can 
now be constructed using 

oo 

V = a n a n H n ((), 

n=0 

where ao = 1 to avoid the coordinate singularity on the shock front, and the remaining 
coefficients arbitrary. In this case, it can be shown that 

2 oo co i 

M = - P - log a - p ]T a n a n H n (() - ]T — a n K n {() 

1 n=0 n=l ZTl 

where 

n-l 

K n{() = a k a n-k [k{n - k)H k H n _ k - (C 2 - 1)-Ff fc _iiy n _jt_i . 

fc=0 

This solution includes additional arbitrary gravitational wave components beyond the 
shock front. However, we note that the curvature singularity which occurs when a = 
is now located on the spacelike hypersurface u + v = in the soliton region and on 
the timelike hypersurface v — u = in the extended region. Since this is a curvature 
singularity, it must form a boundary to the extended space-time as illustrated in figure 6. 
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6 Discussion 



In the above sections, we have considered the physical interpretation and possible exten- 
sions for the soliton solutions with real poles in the case when the seed solution is taken to 
be the Bianchi I vacuum Kasner solution. Very different solutions can be constructed us- 
ing alternative seed solutions. However, the character of the shock front and the possible 
extensions across it is likely to have some similar properties in all cases. 

At least for solitons with a Kasner seed, we have clarified the character of the singu- 
larity that occurs on the shock front, and we have demonstrated a number of possible 
extension across it. The occurrence of thin sheets of null matter, and possible soliton 
extensions, have been discussed in previous literature for some cases. We have clarified 
here the structure of a plane wave extension and given a new explicit vacuum extension 
in which the metric is diagonal. There are clear problems in generalising the non-soliton 
extension to the non-diagonal the equations are then non-linear and superposition 

does not apply. Thus, the full class of permissible extensions has still not been determined. 
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Appendix 

Here we consider whether it is possible to construct an exact one-soliton solution with 
the same global structure as that indicated in Figure 1, but without the presence of thin 
sheets of null matter. Taking the regions I a and h to be the same, the extension would not 
be Kasner. Taking a = \{v 2 — u 2 ) in the one-soliton region I a , it is appropriate to consider 
whether it is possible to choose a gauge such that a = |(w 2 + v 2 ) in the extended region. 
This is clearly C 1 across u = 0, so the Ricci tensor vanishes on this null hypersurface. 

Essentially, we can now prove that solutions of the type outlined above do not exist. 
The proof of these statements are roughly as follows: We start with the general metric 
for a space-time with two spacelike hypersurface orthogonal Killing vectors in the form 

ds 2 = 2e- M dudv - (/(«) + g(v)) ( X dy 2 + X ~\dx - uodyf) 

where the coefficients depend on u and v (or / and g) only. The vacuum field equations 
imply that 



e- M = 4C= e " 5 > (8) 
Vf + 9 

S, = + S , = -l(f + g) ^l±^l. (9) 

X A, 

We may now adopt a gauge such that 

f(u) = \eu 2 , g(v) = \v 2 . 

There is no loss of generality in adopting this form for g(v) and the freedom u — > u'(u) 
has been used to obtain the simplest expression for which /(0) = 0, and f'(0) = to 
avoid any null matter on« = 0. The freedom in rescaling u can be further used to set 
e = ±1. In this case we have 

/' = eu. 
13 



where 



Thus (8) implies that, for M to be continuous across the front u — (on which f + g > 0), 
S must contain the term \ogu. i.e. near the front u = 0, S must behave as 

S ~ log u + const. + . . . 
~ |log|/| + ... 

Thus 

It can then be seen that the first equation in (9) can only be satisfied near the wavefront 
if / < 0. i.e. it is necessary that 

e = -1. 

Thus, solutions with a = u 2 + v 2 do not exist near u — 0. Further, for a vacuum extension, 
it is always possible to choose a gauge such that a = v 2 — u 2 . 
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